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ABSTRACT 

The CORPUS crack growth prediction model for variable-amplitude loading, as introduced by 
De Konrng, was based on crack closure. It includes a multiple-overload effect and a transition 
from plane strain to plane stress. In the modified CORPUS model an underload affected zone 
(ULZ) is introduced, which is significant for flight-simulation loading in view of the once per 
flight compressive ground load. The ULZ is associated with reversed plastic deformation induced 
by the underloads after crack closure has already occurred. Predictions of the crack growth 
fatigue life are presented for a large variety of flight-simulation test series on 2024-T3 sheet 
specimens m order to reveal the effects of a number of variables: the design stress level, the gust 
spectrum severity, the truncation level (clipping), omission of small cycles, and the ground stress 
level. Tests with different load sequences are also included. The trends of the effects induced by 

the variables are correctly predicted. The quantitative agreement between the predictions and the 
test results is also satisfactory. 


INTRODUCTION 

Fatigue crack growth predictions are needed in the aircraft industiy for design and certification 
purposes. Predictions can indicate the trend of the effects of design variables, but it is meaningful 
only if the predictions are quantitatively reasonably accurate. 

Fatigue loads in service, especially for wing structures, are generally a random 
variable-amplitude loading, rather than constant-amplitude loading. Different types of load 
sequences are known to induce load interaction effects, which can result in significant 
accelerations and retardations of fatigue crack propagation. Interaction effects imply that the crack 
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extension in a load cycle depends on what occurred in the preceding cycles, i.e. on the previous 
load history. There are so-called history effects. 

The modified CORPUS model [1] is a modification of the CORPUS model (Computation 
Of Retarded Propagation Under Spectrum loading), which was proposed by De Koning [2] in 
1981. However, the CORPUS model gives a rather conservative prediction for cases, where the 
most severe negative gust is more compressive than the ground load of the flight simulation load 
history. The conservative prediction is caused by the rarely occurring most severe compressive 
gust, which governs the crack opening level in the plastic zone region created by the rarely 
occurring maximum overload. This underload effect of the severe negative gust is lasting very 
long, because it is coupled to a relatively large plastic zone region. In most cases it implies that 
it will last almost during the full flight-simulation history. The modified CORPUS model is 
"uncoupling" underload effects from plastic zone regions created by the maximum overload. In 
addition to the overload affected zone, which is used in the CORPUS model, an underload 
affected zone (ULZ) is introduced in the model.The ULZ is associated with compressive loads 
and reversed plastic deformation induced by those underloads, also by the frequently occurring 
ground load (once per flight). 

The CORPUS model basically is associated with plastic deformation left in the wake of the crack 
(Elber’ s mechanism). Elber [3] observed that permanent plastic deformation ahead of the crack 
tip, generated by high load excursions, is still present on the crack surface after subsequent crack 
growth has occurred. These deformations can cause the crack surfaces to remain in partial contact 
even under tensile loading. Elber introduced the notion of the crack closure and crack opening 
load, defined as the load required for the entire crack surface to be free of contact. Loads 
exceeding this level will lead to crack extension. 

The modified CORPUS model is, to a large extent, similar to the CORPUS model. The 
model includes consideration of plane stress and plane strain condition for the plastic zone size 
and a correction for the effect of high loads on the crack opening stress. The multiple overload 
interaction is maintained as an important feature, which leads to an increasing S^ . 

The model is briefly described and prediction results for a large variety of flight- 
simulation test series are compared to test results. The comparison between predictions and test 
results is restricted to sheet material of the 2024-T3 alloy. Results on 7075-T6 sheet material are 
also covered in [1]. However, 2024-T3 is a more interesting material to check the validity of 
predictions of a crack growth model. First, 2024-T3 is usually adopted for fatigue critical 
components. Secondly, 2024-T3, in view of its larger ductility, gives larger plastic zones than 
7075-T6. As a consequence, larger interaction effects occur, which implies that it is a more 
difficult material to arrive at accurate predictions. 
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FATIGUE CRACK GROWTH MODELLING 


The model will not be described in full detail. The main features will be indicated For a 
complete description the reader is referred to [1] and [2], The process of crack growth modelling 
is illustrated in the flow diagram in Fig.l. As a consequence of crack closure, crack extension 
m every cycle is determined by the effective stress range. In cycle (i) Aa, is- 

Aa 1 = «AK rf „) 

where: AK^ 4 = C AS offj Vita 

The cycle (i) is defined by a maximum S m „ 4 followed by a minimum S^. The effective stress 
range is determined as, 

^eff4 — ^max4 ' ^min4 ^ a op — S min ; 

AS eff4 = S max4 " % if 

where is the maximum crack opening stress level left from the load history. 

Underload Affected Zone 

The underload affected zone should be associated with compressive loads and reversed plastic 
deformation mduced by underloads after the crack has been closed. The K range to be considered 
1S K op - K min- ft may be expected that severe downward loads will be able to induce reversed 
plasticity. It is assumed that this reversed plasticity in the underload affected zone occurs under 

plane strain conditions. The size of the reversed plastic zone is then approximated by the Irwin 
type equation: 
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As long as the crack tip is between a and ARP = a + D u the underload will be effective. Here 
a is the crack length at the moment that the underload was applied. 

An underload affected zone can overlap with another underload affected zone. The zones should 
be stored m the memory of the computer. However, the most severe underload is considered to 
be the dominant one. It is used to determine the S level 

op w * 

Overload affected zone 

The effect of an overload is effective as long as the crack tip is in its plastic zone. The plastic 
zone size depends on the state of stress at the crack tip, plane strain, plane stress or a transition 
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between those two conditions. De Koning [2] derived a special equation for the plastic zone. It 
was based on the Irwin type equation, but it was modified to account for large zones if the stress 
level approaches the net section yield limit. The plastic zone size is: 
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where b is the semi-width of the specimen and y depends on the state of stress. For plane stress 
y = 1/1.32 and for plane strain y = 1/9. 

In the model, the crack tip state of stress depends on the size of the plastic zone relative 
to the thickness. First, the plastic zone size is calculated under plane stress conditions (D ss ). If 
D ss > 0.5t (t=sheet thickness) then the plastic zone is supposed to be in plane stress. If 
Dgj, < 0.35t, the plastic zone is supposed to be in plane strain, size D sn . During the transition from 
plane strain to plane stress (0.35t < D ss < 0.5t), the plastic zone size is given by: 

D = D„ +2D„[(D / 1 -0.35 )/ 0.15 ] 4 (D -D )/ 1 


The overload will be effective as long as the crack tip is between a and ADP = a+D. Also here 
a is the crack length at which the overload occurred. 


During crack growth under variable amplitude loading various plastic zones will be created. If 
they can afreet the crack opening stress in later cycles they must be stored in the material 
memory (plastic zone history, ADP-values labelled as ADPH with H from history). For example, 
consider a crack length a } with a plastic zone size dpj. A second plastic zone dp 2 will be 
memorized if (a 2 +dp 2 ) is longer than (a 1 4dp 2 ). In this case, a new plastic zone penetrates the 
elastic material (ADP 2 > ADPHj). 

Fortunately, not all plastic zones can affect S^, because their effect on is overruled by other 
plastic zones with a higher (S^ > SH^). A higher overload overrules the previous lower 
overload. The series of overload affected zones is characterized by a decreasing series of 
S max -values (SH^^-values). As a consequence only a limited number of plastic zones must be 
stored in the material memory. 


Selection of Crack Opening Stress in Every Cycle 

Two different crack opening occurrences are defined in the modified CORPUS model, i.e. (1) 
crack opening related to the history stress levels (SH^^ and SH^, associated with plastic zones 
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created previously) and (2) crack opening related to the current stress cycle (S ., S • .). The 

history values are related to overloads (SH max ), which produced primary plastic^zonSTand to 

dominant underloads (SH^) of the underload affected zones. As a consequence of the concept 

of the underload affected zone, the dominant (i.e. the highest one) SH max must be combined with 

the dominant SH^, i.e. the lowest SH^ in order to calculate the relevant S (=SH ) for the 
current cycle. op ^ 

The crack opening level induced by the current cycle is determined from S • and the 
successive S^. This is applied in the next cycle only, and only if it exceedTfoe above 

Interaction of Overloads 


s already mentioned in the introduction, overload interaction effects play an important role in 
the CORPUS model. De Koning recognized that overloads with overlapping primary plastic zones 
will cause an extra increase of the crack opening level, which will give more crack growth 

retardation (multiple overload effect). The crack opening stress level for single overload and 
underload is: 

S op = S max (-0.4R 4 + 0.9R 3 - 0.15R 2 + 0.2R + 0.45) if R > 0, and 
S op = S max ( 0.1R 2 + 0.2R + 0.45) if _ 0 .5 < R < 0 

where R is the stress ratio (S^S^). Newman [4] demonstrated that the opening stress level 
oes not only depend on and S max , but also on the maximum stress in relation to the yield 
stress, i.e. S max /a y . In Newman’s analysis elastic-perfectly-plastic material behaviour was assumed 
with a kind of an average yield stress: <r y = (a o 2 +a u )/2. De Koning [2] defined a correction 
ction h wh ich is a good fit to Newman’s results. The correction function h is: 


h = 1 - 0.2 (1 - R)- 


1.15o. 


The corrected stress opening level then is SH op . If interaction between overloads occur, this 
opening stress level will increase due to the multiple overload effect. The opening is increased 
^ter each overload until it has reached a certain upperbound. The equations proposed by De 
Koning will not be reproduced here, since they need a fairly extensive explanation. It includes 

one material constant 5, which was obtained by analyzing empirical data. A foil description is 
given in [1] and [2]. 
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PREDICTION RESULTS 


Flight-simulation test results used in this chapter are adopted from different references. It includes 
results obtained with six spectra, i.e. CN-235 spectrum [5], F-27 spectrum [6], TWIST [7], 
miniTWIST [8,9], FALSTAFF and MiniFALSTAFF[10], and a simplified flight-simulation 
loading [11]. Most of the spectra are related to civil transport aircraft, except FALSTAFF and 
miniFALSTAFF, which are standardized flight simulations for fighter wing structures. The 
material was 2024-T3 Alclad sheet material, except for the tests with FALSTAFF and 
miniFALSTAFF, where 2024-T3 bare material was used. 

Variables of the test programs were: 

CN-235 : design stress level and truncation level 

F-27 : 12 combinations with different 

- gust spectrum severity 

- ground stress severity 

- design stress level 

TWIST : truncation level 

miniTWIST : ground stress level 

FALSTAFF : design stress level 

miniFALSTAFF : design stress level 

The variations of the gust spectrum severity of the F-27 were obtained by increasing or 
decreasing the ratio between gust amplitudes and the mean stress in flight (Sq/S^). There are 
three gust spectrum severities, viz. severe, normal and light. Also three ground load levels were 
adopted with Sg/S^ vary from -0.5 (severe), -0.234 (normal) to +0.125 (light). 

In the miniTWIST and CN-235 tests, the crack increments in the most severe flight (type A) 
were recorded. 

The load sequences used in the simplified flight-simulation loading tests [1 1] are shown 
in Fig.2. During one test all flights were equal. There are three types of flights. The mean stress 
in flight was 80 MPa and the stress amplitude was 40 MPa. Two values of numbers of cycles 
per flight were used, m=5 and m=100 cycles. The purpose of the tests was to study the effect of 
periodic overloads and underloads. Furthermore, a load sequence effect might occur in view of 
the low-high and high-low sequence of the flight profiles II and III respectively. 

Prediction results compared to test results are presented in Figures 3 to 9. 

DISCUSSION 

For most flight-simulation tests the modified CORPUS model and the CORPUS model gave 
approximately similar predictions, with some noteworthy exceptions. The predictions compared 
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to test re S ults diseii 8 ^ below are those obtained with the modified CORPUS model. Predictions 
of both models are given if significant differences were found. 

Effect of spectrum severity 


The effect of the seventy of the spectrum is illustrated by the results in Fig.3. It clearly shows 
an increasing crack growth life if the spectrum becomes less severe. This applies to each of the 

three ground sness levels. The trend is a logical one, and as shown by Fig.3, die trend is 
accurately indicated by the prediction model. 


Effect of design stress level 


The systematic effect of the design stress level is illustrated by Figs 4 and 6 for three different 
oa spectra. The design stress level is represented by either S m „ of the load spectrum or the 
ean stress m flight, S rf . It should be understood that all stress levels of the load spectrum are 
proportional to or As shown by the results, the trend of the effect of Ihe 

stress level was correcdy predicted, although an accurate life prediction was not always obtaLd. 

Effect of the truncation level 


The of truncation of the rarely occurring veiy high loads of a load spectrum (clipping) has 
been well known for a long time (survey in [12],. Predictions in Fig.5 confum this trend shorter 
crack growth lives for lower truncation levels. This is the reason why high truncation levels 
uld be avoided in full-scale flight-simulation tests in order to obtain conservative results. 

Effect of o mitting small cycles 

fi! a gu^t load spectrum cycles with a small amplitude are quite numerous. If such cycles can be 
omitted from a test, savmg of testmg time is considerable. MiniTWIST was derived from TWIST 
^a drastic reduction of the small cycles. Ihe average number of cycles per flight is 100 for 
rf . m . or mmiTWIST. It ieads indeed to a shorter testing time, but the crack growth 
e m ghts is increased about 2 times. The small cycles of TWIST were still damaging. This 
trend agrees with predictions (results not presented here). 

- CA !" Kg ‘ 6 a C ° mpariSOn is made between FALSTAFF and miniFALSTAFF. In 
m AFF about 50% of the smaller load excursions are removed by a rain-flow procedure 

The predictions indicate a negligible effect on the crack growth life in agreement with the test 
results . 
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Effect of the ground stress level 


The ground stress level (Sg,.) in a wing, as compared to the mean stress in flight (S^), is 
dependent on the aircraft configuration (engines, landing gear, fuel tanks) and the rib station. As 
shown by Figs 3 and 7, a more severe ground stress level can imply a significant reduction of 
the crack growth life. This trend is predicted by the modified CORPUS model, but not always 
by the CORPUS model; see next paragraph. 

Comparison between predictions of CORPUS and the modified CORPUS model 

As pointed out before, the CORPUS model and the modified CORPUS model gave similar results 
in most cases. However, this is not true if a rarely occurring negative gust load is a more severe 
down load than the frequently occurring ground load. This condition occurred in the F-27 tests 
with a "mild" ground stress (Sg/S^ = +0.125); see Fig.3c, and in the miniTWIST tests with 
Sgj/Sjjjf = 0 (Fig.7). The predictions of the CORPUS model are indicated by a black dot. The 
CORPUS model gives a significant underestimation, contrary to the modified CORPUS model. 
Actually, this kind of comparative result was the argument for developing the modified CORPUS 
model as discussed before. 

Effect of load sequence 

Load sequences variations were applied in the simplified flight-simulation tests (Fig. 2). The 
modified CORPUS model predicts a negligible difference between load sequences II and III. That 
is in agreement with the test results, see Fig.8. However, the CORPUS model predicts a 
systematic effect if there are 100 cycles in a flight. The life is significantly shorter for sequence 
TIT with the OL at the end of the flight. The different predictions of the two models are caused 
by different S^ developments due to the underload affected zone in the modified CORPUS 
model. 

Crack increment in the most severe flight 


In some tests the crack increment Aa in the most severe flight (type A) could be determined by 
measuring the crack length at the beginning and the end of the flight. Results for miniTWIST 
and CN-235 shown in Fig.9 indicate that the measured increments were significantly larger than 
predicted for those severe flights. The underestimation of the prediction does not have a 
significant influence on the predicted crack growth life, because the more severe flight is a rare 
occurrence. However, it shows a weakness of the prediction model. This is discussed in more 
detail by De Koning et al.[13] (see also [14]). 


554 



Quantitativ e accuracy of the predictions 


The quantitative accuracy of the mdividual predictions was quite good for a fatigue prediction 

Ts7 TT m aTJ? 8 ' ^ average Iati ° ° f predicted clack S™" 1 * life and test life was 
°.8 7 for die modified CORPUS model (85 test results, including insults for 7075-T6) with a 

standard deviation of 0.182. For CORPUS die average ratio was 0.89, practically die same value 

f :r“™ us ’ with standard “ of °- 275 ' wh * h is **• - 


No n-interaction predictions (Miner approach^ 

The non-interaction prediction model is by far the most simple model, because i, ignores any 
effect of previous load cycles. There is no material memory. In each cycle Aa is calculated with 
AK - CASVlta and AS following from S_ and the dfrecdy following S,„„. It is a kind 6f a 
met approach, because the Minter rule also ignores effects of previous load cycles. For the 
&ght simuletion tests with a gust spectrum the average prediction to test result was 0 21 for 
Specimens (29 ,est res " lK >- 11 Mustiates that large interaction effects did occur For 

“ 6 <2 U ° ,eS ' reSUl,S) ,his ra,io 0-45, a higher value associated with smaUer interaction 

effects m the lower ductility alloy. For the manoeuvre spectra the average ratio was 0.45 (2024- 

results only), which also indicate less effective interaction effects. It is associated with the 

effect that manoeuvre spectra have relatively more high amplitude cycles and less low amplitude 
vycics. 

One trend is extremely poorly predicted by the non-interaction procedure, i.e. the effect 
of the truncation level. According to the non-interaction concept a lower truncation level leads 
to a negligible increase of the crack growth life, whereas in reality it implies a significant life 
reduction as shown by test results and is in agreement with the model predictions. 

CONCLUSIONS 

Test results on fatigue crack growth in 2024-T3 sheet material under flight-simulation loading 
were compared to predictions with the modified CORPUS model in the previous chapter Hie 
mam conclusions are summarized below: 

1. Empirical trends of the effects of test variables are correctly predicted. It includes the 
effects of the gust severity of the spectrum, the design stress level, the truncation of rarely 
occurring high loads (clipping), the omission of small cycles, and the ground stress level. 

t ough the CORPUS model also predicts most trends in a similar way, the effect of 
rarely occurring severe negative gusts combined with a not severe, but frequently 
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occurring ground load is not correctly predicted by the CORPUS model. Modified 
CORPUS gives correct predictions also for this case. 

3. A load sequence effect in simplified flight-simulation tests was predicted by CORPUS, 
whereas modified CORPUS predicts this sequence effect to be absent. The latter 
prediction is in agreement with the test results. 

4. The crack extension in the rarely occurring most severe flight was poorly predicted by 
both the modified CORPUS and the CORPUS model. 

5. The quantitative accuracy of the modified CORPUS model was generally good. 
Non-interaction predictions gave insignificant underpredictions and misleading indications 
on the effects of test variables. 
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Fig.l Flow diagram of the Modified CORPUS Model. 
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Fig.8 Effects of overload stress level, S & and the load sequence on the crack growth 
fatigue life. Comparisons to predictions with the modified CORPUS and the 
CORPUS model (•). 
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Fig.9 Crack growth increments in the most severe flight A. Test results and predictions. 
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